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LONG PERIOD PERTURBATIONS 
OF 

EARTH SATELLITE ORBITS 
by 

K.C . Wang 


1.0 INTRODUCTION 


In reference 1, Scheifele and Graf introduced a complete 
first order solution for the orbital motion of a satellite 
perturbed by earth oblateness. This solution was expressed 
in the DS<{> elements. In reference 2, Bond and Scheifele 
expressed the first order short period and secular J 2 sol- 
ution in the non-singular PS4> elements. This theory was 
implemented in an operational computer program named ASOP 
described in reference 3. In references 4 and 5, the PS<J> 
analytical theory was updated to include the drag effects. 

In reference 6, the theory was developed to account for 

the time dependent gravitational harmonics. The drag and time 

dependent geopotential terms have also been included in ASOP . 

Bond also extended the PS<J) theory to include the first 
order long period terms and second order secular perturbations 
due to J 2 , J 3 , and J,. . However, no documentation of 

the equations was ever published. In reference 7, Mueller 
developed a recursive theory to include the first order long 
period terms and second order secular perturbations due to 
zonal harmonics of any order. Mueller's theory plus the 
second order J 2 theory developed by Bond have now been im- 
plemented in ASOP . 
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The purpose of this report is to document all the equations 
involved in extending the PSijs solution to include the long 
periodic and second order secular effects of the zonal harmonics. 
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2.0 METHOD OP SOLUTION 


2.1 Notation 


The DS<J> elements are a set of eight variables which have 
the following description: 

Angle Elements: 

true anomaly 

= g argument of pericenter 

- h longitude of ascending nodp 


tXj = <J> 


a, 


a . 


a. 


= 9, 


Action Elements 

e 2 = G 
B 3 - H 
B, - L 


time element 

related to two body energy 
total angular momentum 
z-component of the angular momentum 
total energy 


These may be canonically transformed to the PS4> elements 
by the following relations: 


1 


= <t> + g + h 

= - /2 ( 4> - G) a in (g + h) 

= - /2(G - H) sin (h) 

- v. 

= 4 > 

= /2( 'I' - (i) r.os (g + h) 

= «/2(G - H) cos (h) 

= L 


( 1 ) 
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The DS4» Hamiltonian for the zonal oblateness problem is 
given by : 


F = Fq + cFj + e F 2 


( 2 ) 


where 


F » $ - 
o 


(two body contributions) 


2L 


F = — 
1 qr 


ol (J 2 contribution) 


W- 

N . . 

- E £ — - — - p (higher zonal harmonics) 

q 11 “ 1 n n- 1 n\ ri 


F 2 =" 


i = .1 iilt 

<5 2 c 


P n are the Legendre polynomials, is the mean equatorial 
radius of earth, and J 0 and !? are oblateness coefficients. 

t n 


2.2 Solution Algorithm 

Von Zeipel's method of elimination of the short and long 
periodic terms is used. The solution first requires the trans- 
formation to eliminate the short periodic terms due to J 2 . 

The generating function is assumed to be of the form 

S = S + eS, 
o 1 

S o give the identity transformation, S 1 is so chosen that 
the new Hamiltonian is no longer a function of short period 
variable <j>. The Hamiltonian has the form 
F'(B\g') = F’ + eF; + c 2 F 2 

A more thorough discussion of the elimination of short 
periodic terms can be found in reference 2. 

An additional transformation must be made to eliminate the 
long periodic terms from F' . This transformation is defined 
by the generating function 

♦ * * 

S ■ S 0 + eS, 
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* * 

Again, S q gives the identify transformation, S 1 is chosen 

such that the long period variable g’ is eliminated from the 

Hamiltonian. The new Hamiltonian has the form 

F"(0 m ) - FJ + eF'j + e 2 F£ 

A more thorough discussion of the elimination of the J 2 
and higher order zonal perturbation long periodic terms can 
be found in references 1 and 7. 

The solution algorithm can be divided into three steps: 

(1) Initialize the primed variables 


P s , as \ 

a ’ k * 0 = ° k ’° + 'Wo 


>"k.O 

Sf, k.0 

3S 

as: 


(3) 


k.O 


■*° V°k.o 3 o k,0/ 


k - 1,2, 3, 4 


(2) Analytical integration of primed variables 



— 

°i,o 

♦ V 


"i 

= 

°2 ,0 

coa( A 2 t ) 

P ' 2 , 0 


= 

°3,0 

<208 ( A jT ) 

P ' 3 , 0 

«; 

s 

°4,0 

+ A. T 
4 



= 

P l,0 



P 2 


p 2 ,0 

cos (A 2 t ) 

°2 , 0 

’J 

- 

P 3,0 

cos ( ) 

+ 0 3,O 



P 4,0 



definitions of A } 

, a 2 , A 3 


stn(A 2 T) 
tiin( A 3 t ) 


oin(A 2 i) 
a in (A^t) 


(4) 


of this report. The relation between time t and the new inde- 

2 

pendent variable x is given by ^ = r /q , the definition of q 
is also given in section 3.0 
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(3) Back transformation 


a 


k 



p k - <>k + 



k - 1 ,2,3,4 


(5) 



3.0 EQUATIONS FOR ELIMINATION OF LONG PERIODIC TERMS AND 
ANALYTICAL INTEGRATION OF PRIMED VARIABLES 

A detailed description of generating function S^ and deri- 
vatives of Sj with respect to the PS$ elements can be found 

in Appendix F of Reference 3. In this section a detailed des- 

* * 

cription of generating function Sj and derivatives of Sj with 
respect to the PS4> elements will be given. The derivatives 
of F” with respect to the DS4> elements will also be given. 


3.1 Generating Function S x 
From Reference 1 we have: 

s i -/WT~ [« - H§ < 2 * 3b ♦ 

\To) q 


6qB)e b sin( 2g) 


■] 


(C) 


S ure terms related to higher order zonal perturbations. 

A 

A detailed description of S can be found in reference 7. 

Now we introduce ,rin( g) and aos( g) 


8in(e) = ^ <o A o, - c,o,) 


6~3 


7 2 


oo 8 (g) = ^ (o 6 o ? + o 2 c 3 ) 
where 

C * ✓2( 4> - G) 

D = S2Tg~ f) 

♦ 

To write in terms of PS elements we introduce the 
following abbreviations 


(7) 


Q =1 


P = 


'P 

1 


J. 

2 


L 


2u 1 t 2 . 2, 

~±ZZ ~ 2 ( S + 

/2oT 2 2 6 


8 


- + ' 


2o 


✓2o, 


( 8 ) 


(9) 


“ ( 1 - p) = 


( 10 ) 
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X - eb^ntn(K) “ ^ q /2<G + H) (a^o^- °7°2^ 
tj> * cb^coD (g) = | 9 /2(G ♦ H) (o 6 o 7 + o 2 o 3 ) 
G - e 2 boin(2ff) * 2x«l* 


Now we have 



(M 


11 

48 


(2 - 3b + 


6qB)©J 


where 

q 



<°6 + 




it 




Ji = |f[f(Sqd + p)(|-b) *b] 

and 

d = (p\i) h 


Let 

T 

a 

T b 

T 

c 

T 


then 


S. 



2 48 


* (2 - 3b + 6qB)T a 
= QT 

b 

A 

= S + T 


T 



(ID 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 


( 22 ) 


Uses****®* 


^hotw*** 
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3.2 Derivatives of S, 


* as 


S lk * 3o. ' k * 1,2,3. ... ,8 
k 

From now on the subscript k represents partial derivatives 
with respect to the 8 PS<f> elements, unless otherwise specified. 


(23) 


# S* T VG J k J 

S lk /9PA T /3F,\ q 


lk HJLffl q J > 

GH k - r[ r * »> <1 - b > + B ] + iK - f (¥ + p ) b . 


+ f(§ - b)p k + t(| - b)2 (dq k ♦ qd £ )| 


(24) 


\ * (H - 1 T ck + (i - 3 v) b k ■ c k + (If + 2 v) * 


bl r k 


ff + 2T b*) *k - I T c - 6BT a°k + 6T » qSB k (25) 


P k = 0 for k = 1 ,3 ,4 ,5 ,7 


‘>2 = 


•2 A* <> 

V 2 


(26) 


P 6 * 


i-ux: 0 
U 6 


— -' A . ■ . 

h) 


for k = 1 , 3 , 4 , 6, 7 


( 27 ) 


2 H j/? 
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G k * 0 for k » 1,3, 4, 7, 8 
G 2 “ " °2 


g 5 - 1 


G 6 * " °6 


(28) 







for k - 1,2, 4, 5, 6, 8 



(29) 


(30) 

(31) 

(32) 

(33) 


# e k = 

2 


" 2o 8 Mp k 1<>r k = 1 ,2,3 ,4 ,5,6,7 


(34) 


- 2o 8 pp - 2p|> 


,/2(G + H) 
°6 <} l ' a 7°2> 


k = ,4,5,8 


Q k“ G G k + 2<G + JO (G k + H k J 


(35) 


X, * (<' 


- i) 


- 0 7 a 2 ) 

/:>( c ' Ti 


/2(G + H) 


2G 


Q* - 


C, G 2 + H) (C 2 * H 2 ) 


1 2G 


J2(G + H) 


(o 6 o 3 - a 7 o 2 y 


2G | Q J 


[ Q J ' G °3 + 2(G Q + H) (G 3 + H 3 ) ] 


+ n 


Q^2(G ♦ ID 


6 2G 
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(a 6 a 3 - Vi) ^ ,8 V 5^(G 6 ♦ y 


oQ^(GTH) 
+ - J -5g 


✓2(G + H) 


X 7 “ ^ a 6°3 " °7 Cr 2^ ~~2G” 


o 0/270 + H) 
- __2 

2G 


7 /2(G + H) 

• *k = ( °6°7 + °2°3 ) 2G 


j^7 ” G G 7 + 2 (G < Vh)* G 7 + H 7 ) ] 


Q k " G G k * 2(G Q + "H) (G k + G k ) ] 


k = 1 ,4, 5,8 


*2 = (0 6°7 + °2 a 3 


(36) 


) S_2 ( G + L 9 g + 9 (Q + H ) 

' 2G [ V 2 G U 2 2(G + H/ U 2 a 2[ 


a 3 Q/2(G + H) 
2G 


*3 = <°6°7 + °2°3 ) 


n S r + Q 
Q 3 “ G °3 2(G + H) VVJ 3 


(G, + H 3 ) 


o 2 Q/2(G + II) 
2G 


*6 = ' V ? + ° 2 ° 3 ) ~ [ Q 6 * 


S 


G °6 + 2( G \ h) (G 6 + »6>J 


£ 


(G* + HJ| 


a 7 Q/2(G + H) 
2G *“ 


*7 = ( g 6 g 7 tC 2 g 3 ) -^ ' ^G ' ^ ' 


G, + 


0 


_ s 

7 G u 7 T 2(G + H) 


(G 


7 + H 7>] 


o,Q/2(G + H) 

O 


2G 
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a 2 

The partial derivatives of S with respect to p , b , e , 
ip and x can be found in Reference 3. 


3.3 Derivative of F ^? with Respect to DS<$> Elements 
From Reference 1, one can find that 


F 2 • 555 6 + 

< - £ (- 3b 2 ♦ 24b - 8) ♦ 18 |£ - i d(j£ ♦ i) 


(37) 

(38) 


o 

(60b 2 - 96b + 32) - Ub(24c 4- 36) 

H is the Hamiltonian of higher harmonics, see Reference 7 
for detailed description. Because the new Hamiltonian is a 
function of only action DS$ elements, from now on the sub- 
script k represents partial derivative with respect to 
those DS<f> action elements. 


2 

6 k = \ [ (e k q “ e2q k )( “ 31)2 + 24b ‘ 8) ] + [(“6b 4- 24)bJ 

bq, 


q 

36b 


q 

L. 


+ ) 
V 


( V - -r* - ^[<Vr + + d <^- p '■ e2p ^ <39) 

P 

» (60b - 96b 4- 32) 4- d(~ 4- ~!) ( 120bb - 96b,)} 

P M k k'l 


- (24c + 36) (B h + Bb ) - 24o 2 Bb 

k k k 


where 

B = 0 

n - 6H 2 

B 2 " " ,4 


* „ 

3 G 3 


( 40 ) 


B ■ 0 
4 
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d 3 “ 0 * 

d 4 - - w<2L)‘ 3/2 

Pi " ' 

P 2 “ - ?! 

P 3 = 0. 

Pi, - -2(pu) !, (2L)' 3/2 

2 2L 

e l " ' H Pi 
2 2 

e 2 = 

e 3 * °- 
2 2 

e 4 = - ~(p + Lp 4 ) 


q L = -0.5 

q 2 = 1.0 
q 3 - 0 

q 4 = - 0.5y(2L)“ 3/2 




(41) 


(42) 


(43) 


(44) 


(45) 


( 46 ) 
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F 2k 


258 < 2 V + «*> + i <k - r 8) 


ft. 


3& _ + 

3F P k + 


aft 

3e 


P 2 aft h 

2 e k 3b b k 


(4? 

(48) 


where 


f 2 i 

f l = \T ( 2 MP + 2q ^ up) 

f 2 

f 2 = ~ U~~( Mp + 2q * /pp) 


f l " ° 

r. = 


\ 2 i 

, = - — rry (vj MP + 2<| /mp) 
4 ( 2 V 


(49) 


Now the abbreviations A t> A^ , , A^ in the expressions 

of analytical integration will be given: 

2 

A 4 = | f 4 (b - 2/3) + u( 2L)‘ 3/2 + —g f (2f 4 6 + f« 4 ) 


E 2 ,A q 4 A 
+ — (H - — H) 
q 4 q 


(50) 


A i - f fb 3 + 255 f2 s + r lft 3 * r 


(51) 


- I [M b - 2 / 3 ) + * b ,] + 2§8 f < 2f ,« + «, > 


+ fr (1} 2 - r + a 3 


(52) 


a. = i + ~ f,(b - 2/3) + ( 2 f x 6 + r,5 l ) 


2 1 

2 * A. 

q 1 “ q~ ' "2 


+ ?- (S, - — H) + A 


( 53 ) 



4.0 CONCLUSIONS 


The equations described in this report have been imple- 
mented into the ASOP program. The program has been checked 
out and verified with results documented in reference 8. Com- 
parisons with numerical integrations show the long period 
theory to be accurate to within several meters after 800 
revolutions. The extension of ASOP to include the long 
period terms, allows the solution to maintain a high degree 
of accuracy even for extremely long prediction intervals. 


ffKCEDWQ 
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APPENDIX 

COMPUTATIONAL PROCEDURE 

The computational procedure for elimination of long peri- 
odic terms and analytical integration of primed variables are 
(^escribed below. First subroutine LONGPP(NN) (long period 
perturbations) is called with parameter 0 , it will return 
initialized primed variable. During the procedure subroutine 
DETERM is called to compute terms related to the higher order 
harmonics. Then subroutine will be called again with para- 
meter 1 , this time it will return the partial derivatives 
q 1 primed Hamiltonian with respect to the DS<J> elements. 
During the procedure subroutine FPRIME is called to compute 
derivatives of higher order harmonics. The sequence of com- 
putation will be given below. The left column gives the 
quantity to be computed, and the right column references the 
equation number in the text. 

LONGPP(O) 

From Equation 

( 20 ) 

( 18 ) 

( 10 ) 

( 19 ) 

( 12 ) 

( 13 ) 

( 14 ) 

subroutine DETERM 

( 30 ) 

( 31 ) 

( 32 ) 
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Computating Sequence (continued) 



o'(0), 


P'(0) 


d 

B 


e 



LONGPP(l) 


From Equation 
(continued) 

( 33 ) 

( 34 ) 

( 24 ) 

( 35 ) 

( 36 ) 

( 25 ) 

( 23 ) 

( 3 ) 


( 20 ) 

( 18 ) 

( 10 ) 

( 40 ) 

( 41 ) 

( 42 ) 

( 43 ) 

( 44 ) 


H, 


A, , A, 


F 2k 

> ,A 3 ,A A 


subroutine FPRIME 
( 39 ) 

( 47 ) 

( 50 ) - ( 53 ) 



